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Abstract 

We consider the planar three body problem of planetary type 
and we study the generation and continuation of periodic orbits and 
mainly of asymmetric periodic orbits. Asymmetric orbits exist in the 
restricted circular three body problem only in particular resonances 
called "asymmetric resonances". However, numerical studies showed 
that in the general three body problem asymmetric orbits may exist 
not only for asymmetric resonances, but for other kinds, too. In this 
work, we show the existence of asymmetric periodic orbits in the el- 
liptic restricted problem. These orbits are continued and clarify the 
origin of many asymmetric periodic orbits in the general problem. 
Also, we illustrate how the families of periodic orbits of the restricted 
circular problem and those of the elliptic one join smoothly and form 
families in the general problem, verifying in this way the scenario de- 
scribed firstly by Bozis and Hadjidemetriou (1976). 

keywords three body problem (TBP), periodic orbits, bifurcations. 

1 Introduction 

The three body problem (TBP), consisting of three bodies with masses mo, 
mi and mi that interact with gravitational forces, has been widely studied 
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in the literature. It is considered either as a simple dynamical system with 
complex behaviour or as a model for explaining the evolution of planets 
and other celestial bodies. In this work, we consider the planar case, where 
all bodies move on the same plane and the indices 0, 1 and 2 refer to the 
particular body. The simplest, yet not trivial, version of TBP is the circular 
restricted TBP (CRTBP), where we assume that one of the bodies is massless 
(say m 2 =0) and the two massive bodies (m ^ 0, m\ ^ 0), called primaries, 
revolve in circular orbits on the plane. When the primaries revolve in elliptic 
orbits of eccentricity e p (=eo = ei) and period T p (=Tq = T\), we have the 
elliptic restricted TBP (ERTBP). 

When the massless body is replaced by a massive one, we get the general 
planar model of the TBP (GTBP). Especially, when one body has very large 
mass in comparison with the other two bodies (say m ^> mi, m Q ^> m 2 ), we 
get the so-called "planetary" GTBP. Continuation and existence of periodic 
orbits in this problem has been studied many years ago by Hadjidemetriou 
(1975) and Bozis and Hadjidemetriou (1976) and recently their results found 
a fruitful field of applicability in the dynamics of resonant extrasolar systems 
(e.g. Psychoyos and Hadjidemetriou, 2005; Haghighipour et ai, 2003; Ferraz- 
Mello et ai, 2005). 

In most cases, periodic orbits are associated with resonant planetary mo- 
tion and their determination in the general TBP can be based on the unper- 
turbed circular problem (i.e. both planets are massless bodies). However, 
not all solutions can be constructed in this way. In Voyatzis et al. (2009), it 
is shown that the main set of families of periodic orbits in the 2:1 resonance 
can be constructed by the continuation of periodic orbits from the restricted 
problem to the general one. In this scheme, the role of families of periodic 
orbits of the restricted TBP is crucial. 

The origin of periodic orbits in the circular restricted TBP can be assigned 
to periodic orbits of the unperturbed problem. The structure of families 
of symmetric periodic orbits is well-known (see e.g. Bruno, 1994; Henon, 
1997). Previous studies (e.g. Beauge, 1994; Voyatzis et al, 2005)) indicate 
that asymmetric periodic orbits exist only in resonances of the form S. = 1, 
q G N, called asymmetric, where for q > 1 the massive body (mi ^ 0) is the 
inner planet and the massless body (m 2 — 0) is the outer one. 

In the elliptic restricted TBP, families of periodic orbits bifurcating from 
periodic orbits of the circular problem also exist having period T = jT v , 
where k and £ are prime integers (Broucke, 1969). Many periodic orbits 
for the elliptic model associated to the dynamics of asteroids and Kuiper 
belt objects have been computed (e.g. Hadjidemetriou, 1999; Voyatzis et ai, 
2005). All the periodic orbits found were symmetric. However, Voyatzis and 
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Kotoulas (2005) found many bifurcation points along the families of periodic 
orbits and conjectured the generation of families of asymmetric orbits. In 
this paper, we compute and show the existence of such families in the elliptic 
restricted model and examine their continuation to the general model. 

In the following section, we present our model and some basic notions on 
the description, existence and continuation of periodic orbits in the particular 
system. In section 3, we show how families of asymmetric periodic orbits can 
appear in the elliptic model and in section 4, how such families are continued 
to the general problem. Finally, in section 5, we conclude our results. 



2 The model of TBP and periodic orbits. 
2.1 The model 

Considering the planetary TBP, consisting of a sun S of mass mo and two 
planets Pi and P 2 of masses m\ <C m and m 2 <C m , respectively, we define 
a rotating frame of reference in the following way (see Hadjidemetriou, 1975) 
: O is the center of mass of S and Pi, the axis Ox is defined by the direction 
S — Pi and the axis Oy is vertical to Ox. Consequently, in the rotating frame, 
the position of the system is given by the coordinates x\ (for Pi) and x 2 , 2/2 
(for P 2 ), while the rotation of the frame is defined by the angular velocity 
9 of the axis Ox with respect to the inertial frame. The Lagrangian of the 
system can be written in the form 



£ = 7^0+^1) 



^(1 - /i) 2 (f 2 + r 2 £ 2 ) + ^ (x\ + y\ + 9\x\ + yl) + 29(x 2 y 2 - x 2 y 2 ))]+V 
m m \ / 



where m — mi +m 2 +m 3 is the total mass of the system, /i — mi/ (m +mi), 
r — Xi/{1 — /i) is the distance of the Sun and Pi and 

Gm^mi Gmim 2 Gm m 2 



r y/((l - fj)r - x 2 ) 2 + yl y/ (fir + x 2 ) 2 + y| 

We remark that 9 is ignorable and, therefore, the angular momentum L = 
dC/d9 is constant. Thus, the system can be reduced to three degrees of 
freedom. By taking the limit m 2 — > in the corresponding equations of 
motion, we obtain the equations of the elliptic restricted TBP. Furthermore, 
by considering the consistent solution r =const, #=const, the equations re- 
duce to the equations of the circular restricted problem. In all cases, the 
system obeys the fundamental symmetry (Henon, 1997; Voyatzis and Had- 
jidemetriou, 2005) 

S : (t,x,y) ->• (~t,x, -y). (1) 
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2.2 Periodic orbits 



An orbit X(£) = (xi(t), X2(t), y2(t), xi(t), x 2 (t), 2/2^)) of the general planar 
TBP is periodic of period T if it satisfies the periodic conditions: 

Xl (T) = ii(0) = 0, 
2Tl(T) =zi(0), 

x 2 (T)=x 2 (0), y 2 (T) = 2/2(0), lJ 

x 2 (T)=x 2 (0), y 2 (T) = y 2 (0). 

A periodic orbit is symmetric when it is invariant under the symmetry E 
and asymmetric otherwise. An asymmetric periodic orbit A is mapped by S 
to another periodic orbit A' which is the mirror image of A. Subsequently, 
all asymmetric periodic orbits appear in pairs (A, A'). 

The first of periodic conditions and the conditon %\ > (or x\ < 0) 
define a Poincare section in phase space. Studying the system on this section 
an asymmetric periodic orbit is assigned to a fixed point and is represented 
by a point in the 5-dimensional space of initial conditions 

IT 5 = {( Xl (0),x 2 {0),y 2 {0),x 2 {0),y 2 {0))}. 

A symmetric periodic orbit crosses perpendicularly the axis Ox i.e., addition- 
ally to the condition ±i(0) = 0, we have y 2 (0) = and x 2 (0) = 0. Therefore, 
a symmetric periodic orbit can be represented by a point in the three dimen- 
sional space of initial conditions 

n 3 = {(x 1 (o),x 2 (o), 2 / 2 (o))}. 

The above periodicity conditions and the spaces H5 and II3 refer either to the 
GTBP or to the ERTBP. Considering the ERTBP, where the semimajor axis 
a\ of Pi is fixed by the normalization of units, the variable X\ is given by the 
eccentricity of the planet Pi, namely Xi(0)=ai(l — e\) or xi(0)=ai(l + ei) 
when Pi is at pericenter or apocenter, respectively. 

In the CRTBP, x\ is constant, defined by the normalization that is 
adopted for the system and x\ = 0. Thus, by considering now the section 
y 2 = 0, a symmetric orbit is defined as a point in the space of initial condi- 
tions n 2 ={(x 2 (0),?/ 2 (0))} and the asymmetric one in the space Il3={(a; 2 (0), 
± 2 (0), j/ 2 (0))}. 

2.3 Continuation 

A periodic orbit of the CRTBP is continued mono-parametrically (e.g. by 
changing the energy integral) forming families (characteristic curves) in the 
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spaces n 2 or IT3 (symmetric or asymmetric, respectively). Along such a 
family the period T varies. In the non-autonomous ERTBP, the periodic 
orbits ought to have period which is an integer multiple of the period of 
primaries T p , i.e. T = kT p . A periodic orbit of the CRTBP can be continued 
to the ERTBP, if its period T' is also equal to T' = kT p or T' = kT p /£. In 
the latter case the generated periodic orbit is assumed to be periodic with 
multiplicity I. 

If we consider a continuation to the space II5 with parameter the x\ 
variable, which, as we have mentioned above, for the ERTBP is related to 
the eccentricity e\ of the primaries, the periodicity conditions are written as 

Y(T; Y(0)) - Y(0) = 0, Y(t) = {x 2 (t), y 2 (t), x 2 (t), y 2 (t)} (3) 

According to the implicit function theorem, the mono-parametric continua- 
tion to n 5 is possible if the Jacobian of system (j3J) has non-zero determinant, 
namely 

|J 4 (T)| = |A 4 -I 4 |^0, A 4 = ^y, (4) 

where I 4 is the 4x4 unit matrix. In case of symmetric periodic orbits, 
the periodicity conditions (J3J) hold for Y(t) = {x 2 (t), y 2 (t)} and Eq. (j3j) is 
rewritten by replacing index 4 by index 2. In general, A 4 has not a unit 
eigenvalue and, subsequently, (j3J) holds and continuation to II5 is generally 
possible. The cases where |J 4 (T)| = correspond to bifurcation points. 

In the ERTBP, the period T in (j4l) is fixed. In the GTBP, T is the time 
span between k successive intersections with the Poincare plane Xi = and 
varies along the family of periodic orbits. Continuation from the ERTBP to 
GTBP is possible when we vary the mass m 2 and keep fixed the variable x\. 
Also, in this case, we obtain again the periodicity and continuation conditions 
(jBl) and (jlj), respectively. 

In the GTBP, the initial conditions of a periodic orbit correspond to 
particular eccentricity values t\ and e 2 , called osculating eccentricities, for 
the planetary orbits. For the representation of families of periodic orbits in 
the space IT5 it is convenient to use the plane e± — e 2 . In this plane, a family 
of asymmetric orbits coincides with its relative family consisting of mirror 
image orbits. Also, the families of the CRTBP lie on the axis e\ or e 2 , since 
in this model one of the planet, P 2 or Pi respectively, has zero eccentricity. 

3 Asymmetric periodic orbits in the ERTBP 

The generation and the numerical continuation of symmetric periodic or- 
bits in the ERTB have been studied in applications concerning the resonant 
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Figure 1: The 1:2 resonant symmetric family E p of the ERTBP and the 

1/2 

biffurcation of the asymmetric family A Ip . 

evolution of minor celestial bodies (Hadjidemetriou 1999; Voyatzis and Ko- 
toulas, 2005). As far as we know, any asymmetric orbits have been given 
in literature for the ERTBP. We propose that such asymmetric periodic or- 
bits can be obtained after two different types of bifurcation discussed in the 
following. 

3.1 Type I bifurcation. 

Along the symmetric families of the ERTBP the linear stability may change. 
At these cases we have critically stable periodic orbits where det | J 4 | = and, 
generally, det | Ja | ^ 0. Therefore, those orbits should be bifurcation points 
in II5 for asymmetric families of periodic orbits. We note that the change in 
linear stability is a quite common fact in the ERTBP (Voyatzis and Kotoulas, 
2005) and, therefore, the families of asymmetric periodic orbits should play 
an important role in its dynamics. 

r 1 /2 

In Fig. HI a 1:2 family of symmetric periodic orbits E p is presented 



6 




which starts from the bifurcation point Be, as an unstable family of period 
2T p . Point Be belongs to an unstable segment of a symmetric family of the 

1 /2 

CRTBP and family E p starts also as unstable and continues parametrically 
by increasing e%. At e\ = 0.5022 the family changes stability type and 
becomes stable. At this point, Be, we have a bifurcation of the asymmetric 
family A 1 /^, which starts as stable. The family E X J 2 , which corresponds to 
the apocenter initial position of the planet Pi, is unstable for all values of e\ 
and, therefore, it does not show a bifurcation of this type. 

A similar structure as above is met in the 1:3 resonance (see Fig. |2J). 
In this case, we have found two symmetric families of ERTBP that show 
bifurcations of type I. The first family E\^ , bifurcates from the 1:3 resonant 

1/3 

symmetric family of the CRTBP, while the second one E 2a starts from the 
circular family of the CRTBP. Both of them start as unstable, but they 
change stability at points B\e and B2E, respectively. At these points, we 
have the bifurcation of the asymmetric families A 1 /^ and A l J^ 2 . 

Bifurcations of type I are also found for the 3:2 resonance. This resonance 
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is neither "asymmetric" nor "external", i.e. the massless planet P 2 revolves 
inside the orbit of the massive planet. In these cases, the CRTBP has no 
asymmetric families. Particularly, the 3:2 resonance has a stable family of 
symmetric periodic orbits which shows a bifurcation point to the ERTBP 
at e\ = 0.4529 {e 2 = 0). Two families of symmetric periodic orbits, E^ 2 
and E a , originate from this point and both of them start as stable. This 
is shown in Fig. [3j Both families change their stability and the bifurcation 
points, BiEi and B 3 e, for asymmetric families are obtained. These families 
are indicated as A 3 / p 2 and A 32 in Fig. We note that the family E 3 J 2 turns 
again to become stable showing a second bifurcation point B 2 e- This point 
does not generate a new family, but it may be assumed as the termination 
point of family A 3 / 2 . 



3.2 Type II bifurcation. 

Considering the asymmetric families of periodic orbits in the CRTBP, which 
exist for resonances in the form oiT 2 /Ti = 1/q, where the primary planet is 
the inner one (Beauge, 1994, Voyatzis et al, 2005), we compute the period T 
of orbits. As we have mentioned the orbits of these families are represented in 
the space n3={(x 2 (0), £2(0), 7/2(0))}. If T = kT p then they consist bifurca- 
tion points and periodic orbits are generated in the ERTBP by continuation 
after varying X\ or, equivalently, t\. Since at the bifurcation x 2 7^ 0, in gen- 
eral, the orbits after continuation should have also x 2 7^ and, subsequently, 
they are asymmetric. 

In Fig. Hk, we show the variation of the period along the 1:2 and 1:3 
resonant families of asymmetric periodic orbits in CRTBP. In particular, we 
present the difference T — To, where T is the period of the orbits along the 
family and T = 2T p or T = 3T p for the 1:2 and 1:3 resonant families, re- 
spectively. Bifurcations for families in the ERTBP correspond to the periodic 
orbits where T — To = 0. We can observe the existence of two bifurcation 
points for each family. We indicate the two of them and present the asymmet- 
ric families which are generated from them in Fig. Hb. From each bifurcation 
point two asymmetric families originate, according to whether the planet Pi 
is initially at pericenter or apocenter. Thus, for the resonance 1:2, we get the 

1 /2 1 /2 

asymmetric families A IIp , which is linearly stable, and A IIa , which is linearly 

unstable. Similarly we obtain the 1:3 resonant families A 1 /^ and A l J 3 a which 
are stable and unstable, respectively. 
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(a) 



(b) 



Figure 4: (a) The variation of the difference T — T along the asymmetric 
families of the CRTBP. The biffurcation points correspond at T — T = 0. 
(b) Type II bifurcation and continuation of 1:2 andl:3 resonant asymmetric 
families in the ERTBP. 

4 Continuation from the restricted to the gen- 
eral problem 

Considering a periodic orbit of the restricted problem, circular or elliptic, 
we can perform continuation by giving and increasing the mass of the ini- 
tially massless body. The continuation of symmetric periodic orbits from 
the restricted to the general problem has been studied by Bozis and Had- 
jidemetriou (1976). In the paper of Voyatzis et al. (2009) relative computa- 
tions were performed for the 1:2 resonance. In this study, we have extended 
our computations for the 1:3 and 3:2 resonances and present two schemes of 
continuation. 

4.1 Continuation Scheme I. 

We mentioned in section 2.2 that the periodic conditions of the ERTBP and 
the GTBP are equivalent. Also, the determinants |J(T)| of the periodic 
conditions, for either the symmetric or asymmetric orbits, are continuous 
functions with respect to the planetary masses. Consequently, in the general 
case, where \J(T) ^ 0|, all families of periodic orbits of the ERTBP are 
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Figure 5: (a) The families E\^ 3 and A 1 / 3 -^ (gray curves) of the ERTBP (1712 = 
0) and their continuation to the GRTBP at ni2 = 0.003. (b) The same for 
the families eV 3 and A 1 / 3 - 



continued in the general problem by varying the masses. 

Particularly, the families evolved in bifurcation of type I do not change 
their topological structure after giving non zero mass to the small body. In 
Figs. and |5b, we show the shape modification of the symmetric families 
E\^ 3 and E^ 3 and the asymmetric ones A X J 3 X and A l J 3 2 , when they are con- 
tinued to the general problem. We present them in the plane of eccentricities 
for ni2 = and m>2 = 0.003. 



4.2 Continuation Scheme II 

According to Bozis and Hadjidemetriou (1976) all periodic orbits of the 
CRTPB are continued in the general problem except those having period 
T = kT p , where, k — 1, 2, ... and T p the period of the circular orbit of pri- 
maries. These orbits are exactly the bifurcation orbits of the families in the 
ERTBP. Therefore, at this point, the topology of the families is not the same 
with that described in the scheme I, but a gap is formed, as we show in the 
following. 

We present the continuation of periodic orbits at T = kT p by an example 

1 /3 

shown in Fig. El We consider the asymmetric families of the ERTBP A/ 
and A 1 / 3 , which bifurcate from the asymmetric family A l J 3 and are shown 
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Figure 6: The generation of families of asymmetric orbits in the GTBP (G\ 

1 /3 

and G2) after the combination of the family Aq of the CRTBP with the 
families A X ff p and A 1 /^ of the ERTBP. 
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Figure 7: The generation of families of the GTBP in the 3:2 resonance. Gray 
thin curves correspond to the families of the restricted problem (mi = 0) 
and the wider ones correspond to mi = 0.0001. Asymmetric families follow 
the continuation scheme I while the symmetric ones follow the continuation 
scheme II. 

in Fig. Hb for m-i = 0. We denote by B the periodic orbit in A c which 
has period T = 3T p and gives the bifurcation point for the families A X J 3 and 

1/3 1/3 ^1/3 

A Ia . The point B divides the family A c in two segments denoted by A lc 
and A x 2 3 . When m 2 takes a non zero value the family A X J 3 breaks at point B. 
The two segments A 1 ^ 3 and A 1 ^ 3 separate from each other and join smoothly 
the families A 1 /^ and A l J 3 a , respectively. So, we get two asymmetric families 
for m 2 7^ 0, denoted in Fig. E]by G\ and At the point B a gap is formed. 

4.3 Continuation of 3:2 resonant orbits 

The above continuation schemes are present in the internal resonance 3:2 as 
it is shown in Fig. [7J We remind that in this case the primary planet is P 2 , 
i.e. m 2 7^ and m x = in the restricted problem. At the bifurcation points 
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Bie, B 2 e and B 3 e, we have the continuation scheme I by increasing the mass 
mi. The position and shape of the families of the ERTBP and GRTBP differ 
slightly as long as mi is small. These families explain the origin of the 3/2 
resonant asymmetric families found by Michtchenko et al, (2006) using an 
averaged model of the planetary three body problem. In the 3/2 resonance, 
we do not have asymmetric orbits generated from bifurcation points of type 
II. However, the symmetric family of the circular problem, lying on the axis 
e2 = in Fig. [TJ has an orbit of period T = 2T P (point Be). Subsequently, 
we get the symmetric families denoted by E^ 2 and E Z J 2 which are continued 
in the general problem {m<i ^ 0) according to the scheme II. 

We remark that, when the masses increase enough, the families may in- 
tersect and new bifurcations arise in the framework of the general problem 
changing drastically the topological structure of families. These bifurcations 
are called "collision bifurcations" and are studied in Voyatzis et al. (2009). 

5 Conclusions 

We have studied some new types of bifurcation of periodic orbits in the planar 
three body problem consisting of a heavy body (star) and two bodies of 
small or negligible mass (planets). When one planet has zero mass (massless 
body) but the second one is massive (primary body), we have the simplest 
models of the circular (CRTBP) or the elliptic (ERTBP) restricted three 
body problem. Considering these models in an appropriate rotating frame, 
we can compute resonant periodic orbits and then perform continuation to 
the general problem where both planets are massive. 

It is known that periodic orbits in the ERTBP are generated from peri- 
odic orbits of the CRTBP and are continued parametrically, by varying the 
eccentricity of the primary body, e p . These orbits are generally symmetric. 
In this paper, we have shown the existence of families of asymmetric peri- 
odic orbits in the ERTBP. There are two types of bifurcations of such orbits, 
which we called as type I and II bifurcations. In a type II bifurcation, the 
asymmetric orbits of the ERTBP are generated from an asymmetric periodic 
orbit of the CRTBP. It is known, that such periodic orbits in the CRTBP 
exist only in the resonances of the form 1/q and therefore such families of 
periodic orbits exist for the RTBP only in these resonances. However, in a 
bifurcation of type I, the asymmetric families are generated from symmetric 
periodic orbits of the ERTBP. Actually, these bifurcation points are critical 
orbits with respect to their linear stability and can be met in any resonance. 
We presented such families for the 1:2 and 1:3 resonance and also for the 
resonance 3:2, which is not of the form 1/q. 
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The families of the CRTBP or ERTBP are continued by varying the 
planetary masses. So, they become families of the general planar three body 
problem. Their continuation with respect to the planetary masses can fol- 
low two different schemes. The continuation scheme I holds for the families 
related to type I bifurcation in the ERTBP. Namely, these families continue 
smoothly by varying the mass preserving their topological structure in phase 
space for adequately small planetary masses. In the continuation scheme II, 
families of both the CRTBP and ERTBP models involve in order to generate 
a family in the GTBP. 

Following the generation and the continuation of the above mentioned 
families we complete the net of families of periodic orbits in the planar general 
three body model and we explain their origin. Periodic orbits of the TBP 
correspond to the so-called "exact resonances" in the planetary dynamics 
and are important since an exosolar planetary system can be trapped in 
such orbits after a physical dissipation process. 
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